Long range triplet Josephson effect through a ferromagnetic trilayer 
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We study the Josephson current through a ferromagnetic trilayer, both in the diffusive and clean 
limits. For colinear (parallel or antiparallel) magnetizations in the layers, the Josephson current 
is small due to short range proximity effect in superconductor/ferromagnet structures. For non 
colinear magnetizations, we determine the conditions for the Josephson current to be dominated by 
another contribution originating from long range triplet proximity effect. 



PACS numbers: 74.45.+C, 74.78.Fk, 85.25.Cp. 

The coexistence of superconductivity and ferromag- 
netism is very rare in the bulk systems. However, it 
can be easily achieved in the artificially fabricated su- 
perconductor/ferromagnet (S/F) heterostructures. The 
S/F proximity effect is characterized by the damped os- 
cillatory behavior of the Cooper pair wave function in 
the fcrromagnct. This phenomenon leads to the non- 
monotonous dependence of the critical temperature of 
S/F multilayers on the F layer thickness and the real- 
ization of the Josephson 7r-junctions (for a review see 
Rcfs.[l|0). In the diffusive limit, the proximity effect in F 
metal is rather short-ranged due to the large value of the 
ferromagnetic exchange field. This is related with the in- 
compatibility between singlet superconductivity and fer- 
romagnetism. 

Interestingly, the non-uniform magnetization can in- 
duce the triplet superconducting correlations which arc 
long-ranged (on the same scale as for superconduc- 
tor/normal (N) metal proximity effect)^. It exists sev- 
eral experimental indications on this triplet proximity 
effecl^i^. However, the transition from usual to long 
range triplet proximity effect was never observed in the 
same system. 

In the present work, we investigate the conditions for 
the observation of the Josephson current due to a long 
range triplet component under controllable conditions. 
The non-colinear magnetization may serve as a source of 
the long range triplet component. However, it is not pos- 
sible to have the Josephson current due to the interfer- 
ence of the triplet and singlet components. Two sources 
of the triplet components are needed to observe the long 
range triplet Josephson effect between them. Then, the 
simplest experimental realization of such a situation may 
be the S/F'/F/F"/S system with the magnetic moments 
of the F', F" layers non-colinear with the F intcrlayer (see 
Fig. H]). The optimal condition for the triplet Josephson 
current observation is when the thicknesses and du 
of the layers F', F" are of the order of the coherence 
length ^/ in the ferromagnet. Indeed, for large d^jdij, 
the triplet component is exponentially small due to short 
range proximity effect in the layers F' and F" , while for 
very thin d^jda, it is also small. Then, we predict that 
the magnitude of the Josephson current in the structure 
with F layer thickness much larger than will be compa- 
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FIG. 1: Geometry of S/F'/F/F'VS junction. The arrows indi- 
cate non-colinear orientations of magnetizations in each layer 
with thickness d^, d, dn, respectively {L — dL + d + dn). 



rable to that of an S/N/S junction with the same length. 

The similar phenomenon could be observed in lateral 
Josephson junctions made of a nanostructured ferromag- 
netic film allowing control on its magnetic domain struc- 
ture. Then, the described effect would give a much larger 
critical current than the one predicted in S/F/S junctions 
with in-plane magnetic domain walls.— 

Besides, the triplet Josephson effect provides the pos- 
sibility of the and 7r-junction realization due to the 
different orientations of magnetic moments in the F' and 
F" layers. Such effect was revealed in Sp/I/Sp junctions 
where Sp are magnetic superconductors with helical mag- 
netic order separated by a thin insulating (I) layer—. It 
was also obtained in diffusive F/S multilayers with non 
colinear magnetizations in successive F layers^. In this 
case, the triplet Josephson effect is mediated by the in- 
verse proximity effect in the thin S layers. It would com- 
pete with the reduction of critical temperature and gap 
amplitude, but these were not taken into account. In 
Ref. i, an idealized circuit-theory model for the triplet 
proximity effect in an S/F/I/F/I/F/S junction was pro- 
posed. The spatial range of singlet and triplet proximity 
effect was not considered. Our work is somewhat compli- 
mentary to these ones. The question about the concrete 
realization and optimization of the triplet Josephson ef- 
fect was outside the scope of these approaches while it is 
of primary importance in the present study. 

We also provide an analysis of the triplet Josephson 
current in the ballistic (clean) limit. In this case, the 
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singlet component reveals the non-cxponcntial oscillatory 
decay but nevertheless the decay of the triplet component 
is even weaker and it is again possible to observe the 
crossover between singlet and triplet Josephson effects. 

The needed conditions for the triplet proximity effect 
observation in the Josephson current are rather stringent. 
The considered system, if realized experimentally, could 
provide an excellent opportunity to study the crossover 
between the triplet and singlet Josephson effects with the 
rotation of the magnetic moment of any of the F layers. 

Let us now calculate the supcrcurrent through a 
Josephson junction made of a ferromagnetic trilayer at- 
tached to superconducting leads, according to the geom- 
etry depicted in Fig. [TJ We assume that the layers are in 
good electric contact and that the magnetizations in the 
layers have the same amplitude. The exchange field h 
acting on the spin of the conduction electrons is parallel 
to the magnetizations, with following spatial dependence: 



layers: 



h{sm(j)LX + cos(j)Lz), < y < di, 

^(y) ^ { hz, cIl < y < dL + d, 

h{sm 4>rx -f cos (j)Rz), dL + d < y < L, 



(1) 



where dz,, d, and d^ arc the thicknesses of each layer, and 
L = dL + d + da is the total length of the junction. Here 
we adopt the sames axes for space and spin quantization. 

We first consider the diffusive limit, when the mean 
free path is shorter than the widths of the layers and 
coherence lengths. For simplicity, we also assume that 
the temperature is close to the critical temperature of 
the leads. Then, within the quasiclassical theory of 
superconductivitjii^, the current flowing through the 
junction is 
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7rT^ImTr[F*(y)a,F'(y)a,], 



(2) 



U1>0 



where the anomalous Green's function F = Fq + F.cr is 
a matrix in spin space and it solves the linearized Usadcl 
equation in the ferromagnet: 

- DF"{y) + 2uoF{y) + ih{y).{a, F{y)} = (3) 

(in units with h = ks = 1). Here, G is the conductance 
of the junction in its normal state, D is the diffusion con- 
stant of the ferromagnet, uj = (2n -I- 1)ttT are the Mat- 
subara frequencies at temperature T, diu. 



{i=x,y,z) 



are the 



Pauli matrices, and the primes denote derivative along 
y-direction. Depairing currents generated by the orbital 
effect have been neglected in eq. ([3]), as usually done for 
ferromagnetic layers with in-plane magnetization^. 

The Usadel equation ([3]) is solved in the central F layer 
in terms of its values at the interfaces with F' and F" 



Fo{y) ± FM = [Fo(di) ± F,{dL)] 
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(4) 



and Fy = 0, as h has no component along y-direction. 
Here, qo = ^2ujjD and q± ~ yj'^iy ± ih)/D. As the 
amplitude of exchange field is much larger than critical 
temperature T^, we may simplify g± ~ (1 ± i)/^f, where 
£,f = yUJJh is the ferromagnet coherence length and 
is much shorter than superconducting coherence length 
^0 ~ \/D/27tTc. The solutions of cq. ([3]) in the other 
layers, as well as their derivative, should match contin- 
uously eq. ^ at each interface. In absence of interface 
barriers with the S leads, they should also take the val- 
ues F{y = 0, L) = F^-^, where F^-^ = (A/a;)eT*x/2 are 
bulk solutions in the leads. Here, A is the modulus of 
the superconducting gap and x is the phase difference 
maintained between the leads. Close to Tc, the gap van- 
ishes as A(T) = [(87rV7C(3))A:2Tc(Tc - T)]^/^. Here, we 
neglect selfconsistency for the gap equation in the leads, 
as usually done assuming that the width of S electrodes 
is much larger than that of F layers, or that the Fermi 
velocity in F layers is smaller—. 

To proceed further with tractable formulas, wc assume 
that F' and F" layers are thin: dL,dR ^ Then, the 
solution in F' layer varies only slightly with y and can be 
put in approximate form: 

F(y) ~ F{dL) + iy~dL)F'{dL) 

{y~dL? 



dl 



[F{dL)~dLF'{dL)-F% (5) 



which satisfies the boundary conditions at y — and 
y = dL- In addition, it should also solve the Usadcl 
equation. Inserting eq. ([5]) into ([3]), we get: 

^[F(di) - dLF'idL) - F^] + -h.{a, F^} 0, (6) 

where a term loF^ was neglected (as h ^ T). Equation 
([6]) yields the results: 
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(7a) 
(7b) 
(7c) 



provided that dL\F'{dL)\ <C |-F(di)|, as can be checked 
consistently from eq. (|4]) when dL £,f- 

Similar results can be obtained for F{y = dL + d). We 
can now evaluate eq. say at ?/ = di, and wc find 
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FIG. 2; Critical current induced by long range triplet 
proximity effect in S/F'/F/F"/S junction, in units of 
(7rGA(T)2/4erc), for varying length of F' and F" layers, at 
dL = da <^ d <^ ^0, and for different orientations of the 
magnetization in the layers. 



I ^ Ic sin X, where the critical current is: 
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The first term in eq. comes from short range singlet 
(i^o) and triplet (Fz") components of anomalous function 
F . It equals the critical current of an S/F/S junction 
with length Its sign oscillates with varying ratio dj^f. 
In particular, when d^ ^f, 
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Thus, its amplitude is also exponentially suppressed. 

The second term in eq. ([8]) comes from long range 
triplet component (F^) and yields: 

let = -Ic,i{d\d\/Sff) sin0L sin(/>_R, (10) 

where Icn is the critical current in S/N/S junctioni^ 

q^d 



e ^ — ' ShOnM ^ 
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(11) 



In particular, in junctions with length d <^ S^q: Icn = 
(7rGAV4erc). The small prefactor {dldy^j) in eq. (fTUl) 
comes from the simplifying assumption dL,dfi ^ ^/ 
that we used in the calculation. As explained in Intro- 
duction, let would be reduced by the exponential fac- 
tor e^^'^'^+'^^-^/^fa.t dL,dR » Thus, at optimal size 
dL,dR ~ the second term let ^ — Icn sin ipL simp 
is much larger than the first one, I^f, provided that the 
magnetic layers have non colinear orientations. For arbi- 
trary lengths dLidfi ~ the critical current originating 
from long range triplet correlation only, at <C (i <C ^Oi 
was also obtained from eqs. ([2]) and ([3]), see FigO We see 
that the triplet contribution to the critical current may be 
observed on the experiment only in the rather small inter- 
val of the F', F" layers thickness: d/,, d/j ~ (0.5 — 2.5)^/. 

The dependence of the critical current with the orien- 
tations of the magnetizations in successive F layers sim- 
ilar to eq. ([5]) was obtained in Refs. 3^. We note that 



the sign of the long range component of critical current 
can be tuned with these orientations. This component is 
absent in the case of only two layers with opposite^^, or 
even non-colinear magnetizations^^. 

The Usadel equations would easily allow generalizing 
the result ([5]) obtained here, (i) Qualitatively, the above 
result should not rely on the assumption that the tem- 
perature is close to Tc and it would be preserved even 
at smaller temperature, (ii) Barrier interfaces between 
the layers and the leads would decrease both short range 
and long range contributions to critical current.— (iii) 
Equation ^ may also describe the case of a ferromag- 
net with magnetic domains and thin domain walls (few 
atomic lengths). If the domain walls are large, the long 
range triplet contribution will be decreased by the factor 
^f/Sw ^ 1, where is the domain wall width, in anal- 
ogy with the theory of enhanced critical temperature in 
S/F bilayers due to domain- wall superconductivity^^. 

Note an interesting possibility to separate the triplet 
and singlet Josephson effects even for relatively thin cen- 
tral F layer d ~ Indeed if its thickness is around the 
first critical value (37r/4)^/, see eq. (O, the temperature 
variation may serve as a fine tuning and provoke the O/tt 
transitioni^ii^. For the S/F'/F/F'VS system, the singlet 
component would vanish at such temperature and only 
the triplet critical current would be observed. 

We consider now the clean limit. Then, the supercur- 
rent flowing through the junction is now given by: 
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where fn(jj) solves the Eilenberger equation in F layer: 
v-VUiy) + ^^U{y) + ^h. \&, Uy)] = 0. (13) 

Here, v = vn is a Fermi velocity, n is a unit vector, G 
is the Sharvin conductance of the ballistic junction in its 
normal state. In addition, the solution of eq. should 
be continuous, and match with the bulk solution in S lead 
where the electrons come from. That is, fn{y = 0) = F^ 
if Uy > 0, fn{y = L) = F^ if Tiy < 0. Again, we neglect 
sclfconsistency for the gap equation in the leads. 

Solving eq. (fT5|) at < y < and rij, > 0, we find for 

fn{y) = /o + f o- that: 

0^/,) = (A/t^)e-''^/2e-2('^±*'')2'/^«, 
^ " " (14) 



/o ± (sin (f>Rfx + cos 

sin^fl/z - cos0_R^ = 



Then, using continuity of / at ?; = d^ and solving eq. (|13p 
ai dL < y < dL + d, we find: 



/o±/. 

fx 



^g-2(^±./.)to-d,)/.„(^^^ ^ .^^^ COS^l), 

~iasin(l)LSd^e-'^'^^y-'^'^'^/'"y (15) 



where a — (A/a;)e ^^f^e '^^'^^hv ^ and we use short no- 
tations = sm(2hdi^lvy\ Cd,^ = cos{2hdL/vy). Similar 
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solution can be found for / at < 0. The supercurrent 
(fT^ is then conveniently evaluated at y = dj^ -\- d/2 and 
we find I = Ic sinx, where: 
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- QSdiSdH COS(/)L COS(^fl^ 

- SdSd^Cd^ COS 01, - Sd^Sd^ sin (/>i sin (/)^] . (16) 

To proceed further, we assume that d^^dji <C ^/ ^ rf, 
where the ferromagnet coherence length = v/h in 
clean limit is much shorter than superconducting coher- 
ence length ^0 = v/2ttTc. Then, 
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Here, the first term comes from short range proximity ef- 
fect. It coincides with the critical current of clean S/F/S 
junction with length d. In particular, at ^/ ^ d ^ ^o, it 
yieldsi^: 
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The second term comes from long range triplet proximity 
effect and yields (for dL ^ d^ <^ ^ d ^ 
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(19) 

It is small under assumption c?/,, dfj <^ . On the other 
hand, at d^^da ^ the critical current p9|) would 
be suppressed by the factor £,j/dLdR <C 1, due to short 
range proximity effect in F' and F" layers. Again, we 
expect a maximum of critical current at dz, ~ da ~ 
with amplitude let oc — I cn sin (j)^ sin (jjR, where Icn = 
{tt AP' G / AeTc) is the critical current of a clean S/N/S 
junction with d <C ^o- The dependence of the critical 



current on the orientations of the magnetizations in F 
layers is similar to the diffusive case. 

The Josephson current through a half-metal (HM) with 
one spin band only is expected to vanish^»i^. However, 
spin-flip processes taking place at S/F interfaces were 
suggested to promote triplet correlation and induce a fi- 
nite supercurrent through the devic a^^'^^i^° . The qua- 
siclassical theory presented in this work assumes that 
ferromagnetic exchange field is much smaller than the 
Fermi energy. Therefore, it is not well suited to address 
quantitatively the case of HMs, when they are compa- 
rable. Qualitatively, the non colinear layers F' and F" 
with thicknesses of the atomic scale would play the role 
of spin flip scatterers with inverse scattering time r^f 
proportional to spin band splitting h. Then, the order 
of magnitude for the triplet induced supercurrent can 
be obtained from eq. by noting that the reduction 
factor dirfij/^j (up to the log term) is proportional to 
l/(rsfi!^_F)^, where Ep is Fermi energy. It is thus propor- 
tional to the probability for an electron from the minority 
spin band to be transfered through HM by spin-flip pro- 
cesses at the interfaces with the leads. 

In conclusion, we determined the Josephson current 
through a ferromagnetic trilaycr. For colinear (parallel 
or antiparallcl) magnetizations in the layers, the Joseph- 
son current is small due to short range proximity effect 
in supcrconductor/ferromagnet structures. For non co- 
linear magnetizations, we determined the conditions for 
the Josephson current to be dominated by another con- 
tribution originating from long range triplet proximity 
effect. In practice the triplet Josephson current may be 
observed in systems with the lateral layers thickness of 
the order of ^/ only. 

The studied structures offer an interesting possibility 
to study the interplay between Josephson current and 
dynamic precessing of the magnetic moment. Indeed we 
may expect the strong coupling between ferromagnetic 
resonance (or/and spin waves) and Josephson current - 
in particular the additional harmonics generation in ac 
Josephson effect. 

We acknowledge Norman Birge and Louis Jansen for a 
critical reading of the manuscript and useful comments. 
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